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r~] ■ Abstract. We provide a general decomposition of the two- photon state produced 

O-i by spontaneous parametric down-conversion in Hermite-Gaussian modes, in the case 

that the pump beam is described by a Hermite-Gaussian beam of any order. We show 
that the spatial correlations depend explicitly on the order of the pump beam, as well 
^ | as other experimental parameters. We use the decomposition to demonstrate a few 

O"* interesting cases. Our results are applicable to the engineering of two-photon spatial 

entanglement, in particular for non-Gaussian states. 
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1. Introduction 

The spatial correlations of photon pairs have been extensively studied in the context of 
quantum optics and quantum information [TJ. In particular, spatial correlations have 
played an important role in the understanding of fundamental quantum phenomena such 
as non-local interference [21 [3] and entanglement pEJ El El [7J EJ E] . Moreover, these spatial 
correlations may be useful for applications such as quantum imaging [TO ], [TT | [2], [T2 | [T3 | [T^j 
and cryptography [151 HE]- For most investigations, correlated photon pairs are produced 
via Spontaneous Parametric Down Conversion (SPDC) [I], which has served as a robust 
and controllable entanglement source. 

It is possible to explore continuous variable (CV) aspects of spatial correlations, 
and in this case one typically considers the near and far-field variables of some reference 
plane. In this CV context, entanglement and non-locality have been observed using 
second-order Gaussian criteria [U El E ELY] , as well as higher-order criteria that are more 
successful for certain non-Gaussian states [El [18]. It has also been shown that basic CV 
quantum gates can be implemented using spatial variables [19] . CV spatial correlations 
can also be used to produce entangled D-dimensional systems [20l l2T] . 

One can also decompose the fields into a discrete set of transverse modes, such as the 
Laguerre-Gauss (LG) ^2H23H2il2Sll2Sll2Zll2Sll2Sll3D]or Hermite-Gauss (HG) modes 
[3lT f32l 133] IMl 135] [36] . Using the LG modes, the majority of experiments have been 
concerned with entanglement in orbital angular momentum [22], [26] [3T] [33] [38] . In terms 
of the HG modes, entanglement in the mode parity has been exploited for experimental 
violation of Bell's inequality [39]. Propagation of entangled transverse modes through 
atmospheric turbulence [ID] and optical fibers has been studied [H]. Investigations of 
quantum information protocols have relied on the entanglement in transverse modes 
[IH H31 H3]. To perform projective measurements in transverse modes, phase masks or 
plane-wave holograms, along with single mode fibers, are typically used. For the HG 
modes, the phase masks introduce a phase shift according to &rg[H n (x)H m (y)], where 
H n is the Hermite polynomial of order n [35]. Images of the first-order plane wave 
holograms are shown in Ref. [13]. Both of these devices transform a higher-order HG 
mode into the fundamental Gaussian mode, which can be selected using a single-mode 
optical fiber. 

An expansion of the two-photon state into a basis of orthogonal transverse modes 
can be written as 



where \jk; a) is single-photon in some transverse mode described by discrete indices 
j, k. Here o~i and o~i describe the widths of the transverse modes. The orthogonality 
condition is 



Experimentally, one is free to choose the values of 0\ and o-i by appropriate coupling 
into an optical system, such as a single mode fiber [221 H21 EE]- An experimental study 




(1) 



j,k,s,t 



{ab;a\a'b';a} = 5 a ^8 b y. 



(2) 
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of this point has been reported in Ref. [30J. In the case of a Gaussian pump beam, 
under certain conditions this expansion has been studied theoretically for the Laguerre 
Gaussian modes (231 ESI E31 EE1 EH] and Hermite- Gaussian modes [l5j |35j [36]. Refs. 
[3T| 134"] considered higher-order Hermite- Gaussian pump beams under the condition 
that o~x = 02 = V2w p , where w p is the radius of the Hermite- Gaussian pump beam. 

Here we derive the expansion of the two-photon state produced by SPDC for 
Hermite-Gaussian pump beams of arbitrary order, using the Gaussian approximation 
of the phase matching function [4"oT W7\ I2"§] . In sections [2] and [3] we review the two- 
photon state and the expansion in HG modes. Explicit calculations are provided in 
the appendix. In section 0] we discuss our results and some interesting examples. We 
show that, under certain conditions, the two-photon state is similar to a classical optical 
beam, and/or photon number states at the output of a beam splitter. Our results are 
applicable to the engineering of two-photon non-Gaussian quantum states [8| IT8]. 



2. Two-photon state 

Here we will be concerned with spatial correlations only. We thus assume that the pump 
laser and down-converted fields are monochromatic, paraxial and have well-defined 
polarizations. These assumptions are justified experimentally by the crystal type and 
geometry, and by using narrow-band interference filters and polarizers, if necessary. For 
a sufficiently weak cw laser, the two-photon quantum state generated by SPDC is [U [H] 



\^) = JJ rfc li rfc l2 $(<Ji, <te) |qi) |qa> • (3) 

The ket describes a single-photon in a plane wave mode. The vector q = (q x , q y ) 
is the transverse component of the wave vector k. The normalized angular spectrum of 
the two-photon state 5>(qi, q2) is given by 

q 2 ) = V( qi + q 2 ) sine ( £|qi 4 ~ qa|2 ) , (4) 

where V(q) is the normalized angular spectrum of the pump beam, L is the length of 
the nonlinear crystal and K is the wave number of the pump field. 

We now consider that the non-linear crystal is pumped with a Hermite-Gaussian 
beam HG nm , generating a two-photon state \ip nm ). We will denote V nm (q) as the 
normalized two-dimensional angular spectrum of the HG mode. The HG modes are 
separable in x and y coordinates, so that V nm (q x ,q y ) = v n (q x )v m (q y ) , where 

(w^q 2 \ 
2/ ® 

and 



D n = -i n * ———=—. (6) 
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Here H n (x) is the n th -order Hermite polynomial. We have included the width parameter 
w explicitly in the definition of the one-dimensional HG modes. Also interesting are the 
diagonal HG modes [IE], which can be defined as 



n+m 



v n (q + ; w)v m (q~; w) = B(n, m, k)v n+m ~ k (q x ; w)v k (q y ; w), (7) 



k=0 



where q± = q x ± q y and 



B(n, m, k) = \l (7l + ™ , fc) | fc! TTTT K 1 " + *) m ] L n- (8) 
v ' ' i y 2( n+m )n!m! k\ dt k 1 v w 

3. Hermite-Gaussian Mode Decomposition 

Since the HG beams form a complete basis, we can expand the two-photon state (jBJ) as 

oo 

h/w>= £ CjLfliA;;^)^;^), (9) 

j,k,s,t=0 

where we define 

\aP; a) = J dqv a (q x ; a)vp(q y ; a) |q) . (10) 
The expansion coefficients are given by the scalar product 

C jklt = ( st '> a l \ 0**5 Vl\^nm) • (11) 

In Ref. [31], explicit analytical expressions for the coefficients C-j™^ were derived 
for the two-photon state of the form with angular spectrum (1131) . There, it was 
considered that the width of the down-converted modes was o~\ = cr 2 = \/2w p , where 



w 



is the width of the pump field. In Ref. [35], it was shown that for a thin crystal and a 
plane wave pump beam V(q) ~ S(q), the down-converted photons are highly correlated 
in transverse modes, such that Cj° k °J t ~ 5j yS 5 k)t . 

In many experimental conditions, the sine function in ffl3|) can be approximated by 
a Gaussian function, such that 

1 [2L . / L| qi -q 2 | 2 ^ 

sine — « V 00 (qi - q2; 6), (12) 



n \ K V 4 ^ 

where 5 ~ 0.257^L/AK, so that the two functions have the same width at half 
maximum [17]. The validity of the Gaussian approximation has been discussed in Ref. 
[IHJ[29]. In this Gaussian approximation, the two-photon angular spectrum is separable 
in the x and y coordinates, and can be written as 



$(qi, Q2) =v n (qix + q_ix\ w x )v m (q ly + q 2y ; w y )x 
vo{qix - q*x\ S x )v (qiy - q2 y ] 5 y ), 



(13) 
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Likewise, the HG expansion coefficients can be separated as = C^Cffl . In 

Ref. [35], an analytic expression was derived for the case of a Gaussian pump beam 
when the down-converted modes are characterized by widths a = a' = 

Here we consider the general situation in which the down-converted modes are 
described by arbitrary width cr, as well as the case of higher-order Hermite-Gaussian 
pump beams V nm . From these, one can construct the coefficients complete C-^ . In 
particular, for c^ = cr 2 = cr, we find 

TV 
fc=0 

where B(a, b, k) is a coefficient defined in Eq. (JSJ) and 

I(j,a,t,s) = Jdqv t (V2q;j)v* s {q;a). (15) 

In the appendix, we derive an explicit expression for the coefficients Cj^\ Our results 
confirm earlier results that the parity is conserved in SPDC [3U [SI [35], so that 
parity (a + b) = parity(n). Below, we provide some numerical examples and discuss 
some interesting cases. 



4. Interesting Cases 

Jf.,1. Two-mode squeezed state 

By adjusting the width of the collected down-converted modes, one can manipulate the 
terms in the expansion Q). A particularly interesting case is when a = o' = V2w5 and 
the pump beam is a Gaussian. Here the coefficients are given by 

(o) ws (w-sy 

where 8j s is the Kronecker delta. This corresponds to the Schmidt decomposition 
calculated and observed in Ref. |35j : 

oo 

H'oo) = y^,V \ab; cr'} \ab; a') , (17) 

ab=0 

where \f\t = C aa and we specify explicitly that the down-converted modes have widths 
a'. This is similar to the two-mode squeezed vacuum state for intense fields, which 
displays photon number correlations between the modes [50]. Here, the correlation is 
between the HG mode indices in both the x and y directions. For w — 8, only one 
coefficient = = A is nonzero, and the state is separable. This is true only 

in the Gaussian approximation. However, for the wave function given in Eq. ffTB"]) . the 
state displays a minimum amount of entanglement [IHl 
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Figure 1. The coefficients C^j/ for several values of w/S, with a = \j2w5. The 
intensity scale goes from the largest negative (blue) value to the largest positive value 
(red). ' 
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Figure 2. The coefficients C ab for several values of w/S, with a = \/2w5. The 
intensity scale goes from the largest negative (blue) value to the largest positive value 
(red). 
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Figure 3. The coefficients C^/ for several values of w/S, with a = y/2wS. The 
intensity scale goes from the largest negative (blue) value to the largest positive value 
(red). ' 



4-2. Correlated HG states 

For the same width parameter a', and considering now higher-order HG pump beams, 
the two-photon state takes the simplified form \ip) nm = \<fi n ) \4>m), where 

oo n 

i<m = E E c SUr(i«>ii«+'-> 2 + i«+'->ii«> a ). ( 18 ) 

a=0 r=0,par(n) 

and similarly for \<p m ). Here |a) are single photons in one-dimensional HG modes, in 
analogy to Eq. ( ITUj) . For n, m ^ 0, these states are always entangled. Figures (TJ [2] and 
|3]show the coefficients C^) for various values of a' in the cases n = 1, 2, 5, respectively 
Different plots in each figure correspond to different values of w/8, the ratio between 
the width of the pump beam and the phase matching function. Correlation of the form 
given by Eq. ( 1T8|) can be observed. For w/5, the states show the minimal amount of 
entanglement for each value of n, as has been shown in other studies jlEj. In addition, 
the minimal number of HG modes in the expansion increases with n, indicating that 
entanglement also increases with n [31]. We note the alternating phase changes along 
the diagonal of the figures for coefficients of the HG decomposition when w/8 < 1. In 
the case of the Gaussian state described in Eqs. (|T6l) and ffTTj) . the ratio w/S determines 
whether the Gaussian functions of two-photon state presents momentum correlations 
(5 > w, qi ~ q 2 ) or anti-correlations (5 < w, qi ~ — q 2 )- As can be seen in Eq. ( !T6|) . 
when j is odd the sign of the coefficient Cff is negative for correlations and positive 
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Figure 4. a) Left: The probability distribution | <^>3 1 2 given by Eq. (TTS)) . The two 
dimensional coordinate system cooresponds to one spatial dimension q\ and q2 of each 
photon produced by the non-linear crystal (right), b) Analogy of state \<j> n ), given by 
Eq. (jljjp. with a classical Hermite- Gaussian mode with spatial variables q x and q y . c) 
The output state produced when a photon number state |n) is superposed with vacuum 
|0) at a beam splitter is also analogous to two- photon state (JTHJ) . 



for anti-correlations. Thus, for 5 > w the phase oscillates as a function of the sign of 
j, and this is related to the spatial symmetry (correlations or anti-correlations) of the 
two-photon state. For higher order HG pump beams as in Figs. [U - El a similar effect 
is present. In addition, the number of modes present depends explicitly on \w/5\, and 
the coefficient C^ 1 decays exponentially with a and b. Thus, by controlling the ratio 
\w/6\, the number of transverse spatial modes can be manipulated. This leads directly 
to control of the entanglement of the two-photon state jlHl [3H [29] . The states f lT8|) are 
similar to non-Gaussian states studied in Refs. [511 El HE] with n — 1, which display 
interesting properties. We expect states with n > 1 to display additional novel features 
[521. 



For w p = 5 and a = y/2w p 5, the coefficients C t 



ab 



n\/2 n a\b\, and 



« = E 



a=0 



2 n a\(n - a)\ 



\n — a} 1 \a) 2 . 



(19) 



The quantum state describing the spatial variables of the photon pair are analogous to 
other physical situations in classical and quantum optics, as illustrated in Figure HI The 
states \4> n ) in (I19p are then analogous to the usual diagonal HG beams in classical optics 
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defined in Eq. (171) , with index m = 0. In this case, it is straightforward to see that 



TV 

g <"- -*> = ' W-*)r (20) 



giving 



v n (q + ;w)v (q_;w) = 



nl 



k=0 



2 n n\(n - k)\ 
v n _ k (q x ;w)v k (q y ;w). (21) 



The two-photon quantum state ([19]) is thus analogous to the two-dimensional classical 
beam defined in Eq. ( |2TT) . where the one transverse spatial variable of each photon-say 
qi and g2~play the role of the transverse wave vectors q x and q y of the beam. Similar 
two-photon states with this property have been observed for LG modes [5TJ [52] . 

To make an analogy with quantum optics of field modes, the wave function of state 
\<f) n ) is equivalent to combining a photon number state |n) and a vacuum state on a 
50/50 asymmetric beam splitter. The quantum state of the two modes at the output 
ports of the BS is 



,t 



\ \ \ \ Y \ tP dn-j), b'),, (22) 

where is the usual creation operator. Comparison with Eq. f|T9|) shows that the two 
states are equivalent. 

For the general case in which a ^ a', the quantum state still presents parity 
correlations. Figure [5] shows four examples for n — 2. If one chooses a = V2w p 
(graphics on left), in addition to the parity restriction we have a + b > n as can be seen 
by the absence of the a = b = term in the two graphics on the left, in accordance 
with previous work [ST] IM] . Thus, it is possible to produce two-photon states where 
lower order transverse modes are prohibited. The two figures on the right correspond 
to a = y/25. In this case parity correlations can also be observed, but we notice that 
there is also a nonzero coefficient when a = b = 0. Interestingly, when w/5 = 2, only 
the a = b = coefficient is negative. For larger w/5, we observed numerically that 
additional negative coefficients appear. 



5. Conclusion 



We have derived a Hermite-Gaussian expansion of the two-photon state produced by 
spontaneous parametric down- conversion when the pump beam is an arbitrary Hermite- 
Gaussian mode. Use of the Gaussian approximation to the phase matching function has 
allowed us to arrive at analytical expressions for the expansion coefficients. Our results 
show that a number of interesting two-photon states can be produced by changing the 
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Figure 5. The coefficients C^ b for a = \pLw and a = \/28, with w/8 = 2 and 
iu/<5 = 1/2. The intensity scale goes from the largest negative (blue) value to the 
largest positive value (red). 



Hermite- Gauss indices of the pump beam, as well as the widths of the pump beam 
and the phase matching function. An interesting feature of the HG mode expansion is 
the fact that, under certain conditions, one can consider a single spatial dimension for 
each photon. This allows for simplification of the theoretical analysis and experimental 
investigations [35] • We have also shown that in certain cases the two-photon state in 
one spatial dimension is analogous to a diagonal Hermite-Gaussian beam in classical 
optics, and also to a photon number state superposed with vacuum at a beam splitter. 
The analogy between a single spatial dimension and a single harmonic oscillator allows 
for interesting parallels to be drawn. For instance, viewing the spatial correlations 
as forming a "two-photon beam", as illustrated in Fig. H] a), allows one to identify 
interesting features, such as the correlation vortex reported in Refs. [51], [52]. Since 
this image relies on a single spatial dimension of each photon, it is more adequately 
described using the HG modes, described in cartesian coordinate systems. We expect 
our results to be useful in the engineering of spatial correlations, in particular in the 
case where non-Gaussian properties of the quantum state are to be explored [81 fT8l l5T] . 
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Appendix 

The one- dimensional expansion coefficients are Cj£ = (v a \ (vb\ip n ). Using the definition 
of HG modes (jSJ) and the two-photon state (j3J), we have 

v* a (qi; o- 1 )vl(q 2 ; a 2 )dq 1 dq 2 (23) 

Defining 

q ± = -j=( qi ±q 2 ) (24) 

we can rewrite the integral in terms of q± variables. Let us consider that the down- 
converted modes satisfy o\ = o~ 2 = a. Using the definition of the diagonal HG modes 
(J7J), we can rewrite Eq. ( 123]) as 

N f 
C ab = Yl fi ( a ' & ' k ^ / d( l+ v n(\feq+] w p )v* N _ k (q + ; a) x 

k=0 J 

dq_v (V2q_;8)v* k (q_;a), (25) 

where N = a + b and B (a, b, k) is defined in Eq. (J8J). The integrals in Eq. (1251) can be 
solved separately. Then coefficients can then be written as 

N 
fc=0 

where 

J( T ,<7,t,s) = Jdgv t (V2q-i)v* B (q;<r) (27) 

Expanding the Hermite polynomials in power series [53] and integrating each term, one 
has 



, 7^!/°? /2) (-1)- 



7 2 + (T 2 /2 ^ £ , (t _ 2 ^ 

^ j!( S -2j)! X 

• • s-2j 



2^7 V f 2a 



7 2 + a 2 /2/ V7 2 + ^ 2 /2 



x 



^2 2 2 
where T is the usual Gamma function 
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When t = 0, a, 0, s) can be simplified to 

. ^aD D*s\ fa 
7(7, a, 0, s) = ^— r - 1 . 29 
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